Counting statistics of charge pumping in an open system 
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Electron counting statistics of a current pump in an open system has universal form in the weak 
pumping current regime. In the time domain, charge transmission is described by two uncorrelated 
Poisson processes, corresponding to electron transmission in the right and left direction. Overall 
noise is super-poissonian, and can be reduced to poissonian by tuning the amplitude and phase of 
driving signal so that current to noise ratio is maximized. Measuring noise in this regime provides 
a new method for determining charge quantum in an open system without any fitting parameters. 

Electric current through an open electron system, such 
as a quantum dot well coupled to the leads, can be in- 
duced by modulating its area, shape, or other parameters 
0. The possibility to generate a DC current through a 
quantum dot by cycling potentials on the gates was pro- 
posed by Spivak et al. |^ and realized experimentally by 
Switkes et al. Theory of pumping in open systems 
was developed by Brouwer [Q and by Zhou et al. 
Brouwer made an interesting observation that time aver- 
aged pumped current is a purely geometric property of 
the path in the scattering matrix parameter space, insen- 
sitive to path parameterization (also, see Refs. p^ , p^ ). 
Zhou et al. demonstrated |^ that pumping provides a 
new approach to a detailed understanding of mesoscopic 
transport. Recently, a number of issues related to incom- 
plete agreement between theory and experiment were ad- 
dressed P-pUl . An interesting extension of these ideas to 
mesoscopic superconducting systems was proposed p^ . 

In this article we discuss current fluctuations in para- 
metrically driven open systems. In the regime of in- 
terest, called "adiabatic pumping," system parameters 
change slowly compared to transport time through the 
system. This problem is different from adiabatic trans- 
port proposed by Thouless [p^ , involving an open sys- 
tem with a gap in the excitation spectrum. Thouless 
pump is adiabatic in the quantum-mechanical evolution 
sense, provided that the driving frequency / is smaller 
than the energy gap in the system. Current in the Thou- 
less pump is quantized in the units of e/, which has been 
demonstrated in quantum dots [ p^|jr^ and also motivated 
proposals to detect fractional Quantum Hall charge . 
We demonstrate below that, although in an open system 
pumped current is not quantized, charge quantum can 
still be detected from noise measurement. 

Coherent transport through an open mesoscopic sys- 
tem is described ]l9| by a unitary scattering matrix S 
which depends on externally driven parameters, and thus 
varies with time. The matrix S(t), as a function of time, 
defines a path in the space of all scattering matrices. For 
a system with m scattering channels, the matrix space is 
the group [/(to) = SU{m) x U{1). 

In the pumping experiment one can, in principle, real- 
ize any path in the space of scattering matrices. In this 
article we consider the regime of a weak pumping field. 



FIG. 1. Current to noise ratio, I / J = e~^{u ~ v) / {u + v), 
as a function of the driving signal parameters (jisj) for a sin- 
gle channel pump. The right and left transmission rates u 
and V are given by (|2^). Two harmonic signals driving the 
pump are characterized by their relative amplitude and phase, 
w = (^1/1^2)6'". Maximum and minimum, as a function of w, 
are I/J = ±e~^, where e is elementary charge. 



when the path S{t) is a small loop. We show that in 
this case the distribution of charge q transmitted per cy- 
cle is fully characterized by only two parameters, average 
charge flow per cycle, I = f{q), and noise, J — /((q^)). In 
the time domain charge transport is described by two un- 
correlated Poisson processes for independent single elec- 
tron transmission to the right and to the left. The gen- 
erating function for charge distribution over N pumping 
cycles in this case is 



X(A) 



(1) 



Here the rates u and v of right and left transmission per 
cycle are given hy u — v — I / ef , u -\- v — J / f , where 
e is elementary charge. Counting probabilities p„ can be 
found from Fourier decomposition, x(A) = ^s^^^Pn- 

Poisson statistics, identical to that of conventional clas- 
sical shot noise, makes it possible to use pumped current 
noise as a new method of measuring charge quantum. 
However, since right and left current fluctuations (^ are 
independent, one needs a way to separate the two Poisson 
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processes. This can be achieved, as discussed below, by 
varying amplitudes and relative phases of external driv- 
ing signals. Either the right or the left transmission rate, 
u or w, can be nulled. The parameters for which this 
happens give extremum (maximum or minimum) to the 
current-to-noise ratio I / J — see Fig.|l|. Once (|^) is re- 
duced to a single Poisson process, this ratio gives charge 
quantum without any fitting parameters, J/I = e. 

So far, only noise in nearly open systems has been used 
to detect quasiparticle charge. In particular, shot noise 
measurements in Quantum Hall point contacts [ pOpl| ] use 
backscattering current of a conductance plateau. Theo- 
retical discussion of ways to detect fractional charge in 
Quantum Hall systems |22| and in Luttinger liquid [ ^J24| 
also focuses on weak backscattering regime. Based on the 
present analysis, we conjecture that the requirement of 
ballistic transport is not necessary if current is induced 
by weak pumping, rather than by a DC voltage. The 
method discussed below allows to determine charge quan- 
tum in open systems with significant scattering. 

Another result we report concerns general dependence 
of counting statistics on the path in matrix space. Differ- 
ent paths S{t)^ in principle, give rise to different current 
and noise. However, there is a remarkable property of 
invariance with respect to group shifts. Any two paths, 



S{t) and S'{t) = S{t)Sa, 



(2) 



where 5*0 is a time-independent matrix in U{m), give 
rise to the same counting statistics at zero temperature. 
We note that only the right shifts of the form (j^) leave 
counting statistics invariant, whereas the left shifts gen- 
erally change it. One can explain the result (||) qual- 
itatively as follows. The change of scattering matrix, 
S{t) — > S'[t) = S{t)So, is equivalent to replacing states 
in the incoming scattering channels by their superposi- 
tions = 'S'o/jV'^- At zero temperature, however, Fermi 
reservoirs are noiseless and also such are any their su- 
perpositions. Correlation between superposition states 
of noiseless reservoirs is negligible, because all current 
fluctuations arise only during time-dependent scattering. 
Therefore, noise statistics remain unchanged. A simple 
formal proof of the result (|2[) is given below. 



I. GENERAL APPROACH 

The distribution of transmitted charge is characterized 
p5| , p6t by electron counting probabilities p„, usually ac- 
cumulated in one generating function xW = J2^^^^Pn- 
The function xi^) is given by Keldysh partition func- 
tion, describing evolution in the presence of a count- 
ing field A which is an auxiliary gauge field having 
opposite signs on the forward and backward parts of 
Keldysh time contour [^,^. By a gauge transfor- 
mation, the time-dependent scattering matrix becomes 



S\{t) = e'^i"^S{t)e for the forward time direction, 

and Si^{t) = e-*^'"3 5't(t)e'^'"3 for the backward direc- 
tion, where (73 is a diagonal matrix with eigenvalues 1 
and —1 for the right and left channels. Then 



X(A)= dei{t + n{t,t'){fx{t)~t 



f^{t)^Sl^{t)Sx{t). n(t,t') 



2Ti{t -t' + iS) 



(3) 
(4) 



where fi is the density matrix of reservoirs at zero temper- 
ature. The determinant of an infinite matrix (^) requires 
proper understanding and definition p7| , p8[ . 

There are two ways of handling the determinant (|^). 
For periodic S(t), one can use frequency representa- 
tion 1 26 1 in which n is a diagonal operator, n(uj) = 
9{—Lo). The operator S{t) has off-diagonal matrix ele- 
ments S^'^uj with frequency change being a multiple of 
external (pumping) frequency, lj' ~ uj — nO. — 2Tmf. 
In this method the energy axis is divided into intervals 
nf2 < Lu < (n + l)r2, and each interval is treated as 
a separate conduction channel. In doing so it is con- 
venient (and some times necessary) to assign separate 
counting field A to each frequency channel, so that the 
field A may acquire frequency channel index in addition 
to the usual conduction channel dependence given by 
in (^. This procedure brings (^ to the form of a deter- 
minant of a matrix with an infinite number of rows and 
columns. This matrix is then truncated at very high and 
low frequencies, eliminating empty states and the states 
deep in the Fermi sea which do not contribute to noise. 

This method was used in Ref. |^ to study noise in a 
two channel problem described by a 2 x 2 matrix 
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be 



-B 



(5) 



which is unitary for \AY 



+ = l,Aa + Bb 



0. Charge distribution for N pumping cycles is 



X(A)- (l+Pi(e»^-l)+P2(e-'^-l)) 



N 



(6) 



with pi = |a|V(|a|' + 1^1') and P2 = 



6|V(|ap + |6p). 



Alternatively, the determinant in (^) can be analyzed 
in the time domain. This representation is beneficial 
when an orthogonal basis of functions can be found that 
provides a simple enough representation of (||). For deal- 
ing with periodically driven systems, it is convenient to 
close the time contour by imposing antiperiodic bound- 
ary conditions on the interval < t < tq = 2TrN/il, 
where TV is the number of pumping cycles. Upon doing 
this, the reservoir density matrix n acquires the form 



n{t, t') = 



2to sin {tt (t - t' + i6)/ To) 



(7) 



The reason for introducing periodicity (^ is that in this 
problem one deals with probability distributions which 
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are stationary in time, so that x(A) is multiplicative, and 
Inx(A) (X N at N 1. Conveniently, the periodicity 
of (|^) makes multiplicative character of xW an exact 
property, true even for iV ~ 1. 

One encounters a number of interesting cases which 
can be handled in the time domain in the problem of 
noise in voltage-driven systems. An external driving volt- 
age V{t) apphed across the system can be incorporated 
in S* as a time-dependent forward scattering phase. This 
is achieved by a gauge transformation making the scat- 
tering matrix time-dependent: 



with T 



(0) 



S-i^M"^ and Ax{t) = e'i''^A{t)e-'i' 



(8) 



The formula defines a circular path in U{m) of radius 
which depends on the system conductance. Full statis- 
tics have been studied for a large family of paths of the 
form (^). The statistics was found to be binomial for the 
DC voltage case as well as for the AC case [ p7||28| 
with a particular time dependence V(t) obtained from 
the criterion of minimal noise. 

An attempt to adapt these results to the problem of 
pumping noise was made by Andreev and Kamenev [ p9[ . 
For several matrix paths S{t) e SU{2) obtained from 
by exchanging incoming scattering channels, bino- 
mial charge distributions arise, not surprisingly, with 
transmission and reflection probabilities exchanged. An- 
other matrix considered in Ref. |29|, r — —r' = cos fit, 
t = t' = sinfii, is related to (5|) by the transforma- 
tion (Q) with 5*0 = ((73 + (T2)/V2 and parameter values 
a,b = I/V2, A,B ^ 0. The resuh (Eq.8 of Ref. (29|) 
in this case agrees with (^, in full accord with the in- 
variance property (||). However, in the context of the 
pumping problem, the paths S{t) considered in Ref. ||29| ] 
appear to be less relevant than, say, the paths (^) in 
voltage-driven systems. In general, the dependence of 
the scattering matrix S{t) on the parameters externally 
controlled in the pumping experiment is not known. Be- 
cause of that, the results for particular paths S are of 
less interest than the properties that hold for sufficiently 
general families of paths. 



II. CALCULATION 

For a weak pumping field we shall evaluate in the 
time domain by expanding lndet(...) in powers of SS and 
keeping non- vanishing terms of lowest order. In doing 
so, however, we preserve full functional dependence on A 
which gives all moments of counting statistics. We write 
S{t) = e"^^*'5^°^ with antihermitian A{t) representing 
small perturbation, trA^A <C 1. Here S''^*'^ is scattering 
matrix of the system in the absence of pumping. Substi- 
tuting this into (0) one obtains 



Tx{t) ^ f. 



(0) 



JTaW = 5i°]te-^-We^^W5f (9) 



Now, we expand (| 



Inx(A) = IndetQo + tri?- itri?^ + itri?^ - ... (10) 

where Qo = 1 + n{f^^^ - 1) and R = Q^^MTx. At zero 
temperature, from = n \i follows that detQo — 1 and 
R = Sf ^ '^h (e-^-A(*)e^^(*) - 1) 5*^"^ Therefore, 



Inx(A) = tr nil - -ix{nMY + -iY{hM f 
2 3 



(11) 



where M 



1. Note that at this stage 
there is no dependence left on the constant matrix 5"*^°), 
which proves invariance under the group shifts (|^). 

We need to expand (11) in powers of the pumping field, 
which amounts to taking the lowest order terms of the 
expansion in powers of the matrix A[t). One can check 
that the two 0{A) terms arising from the first term on 
the RHS of (|ll|) vanish. The 0{A^) terms arise from the 
first and second term in ([l^) and have the form 

Inx = ^tr (n {AI^ + A\ -2A^xAy))-\lT{hBxf 

(12) 

with B\{t) — A\(i) — A^\{t). At zero temperature, by 
using n} — h, one can bring (|l^) to the form 

i tr [n [Aa, A-a]) + I (tr {n^Bl) - trifiBxf) (13) 

The first term of (|l^) has to be regularized in the 
Schwinger anomaly fashion, by splitting points, t',t" = 
tzL e/2, which gives 

n{t', t")tr [A^xinA^it') - Ax{t")A_x{t')) dt (14) 

Averaging over small e can be achieved either by insert- 
ing in ( p^ ) additional integrals over t' , t" , or simply by 
replacing Ax{t) -> \ {Ax{t) + Ax{t')), etc. After taking 
the limit e — > 0, Eq.(p^) becomes 

■^<j>tr{A^xdtAx- AxdtA,x)dt (15) 

The second term of ( |l3|) can be written as 
1 f f tr {Bx{t) - Bx{t')f 



4(27r) 



{t-t') 



l\2 



-dtdt' 



(16) 



Now, we decompose A ^ aQ + z + z\ so that [173, Og] = 0, 
[ct3,z] = -2z, [cr3,^^] = 2zt. Then 



Ax = e-'-""" Ae'-"' 



Bx = ( e*^ - e"** W = - z 



(17) 
(18) 
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Substituting this into (15) and ( |l^ ) one finds that in 
terms of W{t) these two expressions become 



sin A 



tY{[a3,W]dtW)dt 



and 



(1-cosA) j;j;tr{W{t)--Wit')Y 
2(27r)2 



{t-t'y 



dtdt' 



(19) 



(20) 



Hence In x indeed depends on A asw(e^^-l)+z;(e-'^-l). 

Eq. (|^ is essentially identical to the result obtained by 
Brouwer for average pumped current . The integral in 
( p^ ) is invariant under reparameterization, and thus has 
a purely geometric character determined by the contour 
S{t) in U{m). Eq.(pO[) represents a generalization of the 
expression for noise induced by time-dependent external 
field considered in Refs. [|o|j2^]. 

The parameters u and v in (|^) can be expressed 
through z{t) and z''{t) in a simple way. Let us write z{t) 
as z+(i) + z_(t), where z+{t) and z^{t) contain only pos- 
itive or negative Fourier harmonics, respectively. Then 

/.r (.la.,-,, c,,,l)*.i:.. (21) 



fermions, we argue that it can be applied to interacting 
systems as well. Poisson distribution results from the 
absence of correlations of subsequently transmitted par- 
ticles, which must be the case in any system at small 
pumping current. Using the dependence of Poisson rates 
u, V on the driving signal to maximize I/J allows to 
eliminate one of the two Poisson processes (^ and then 
obtain charge quantum in the standard way as e = J/ 1. 

I am grateful to Dmitry Novikov for useful discussion. 
This work was supported by the MRSEC Program of the 
National Science Foundation under Grant No. DMR 98- 
08941 



UJ>0 



w = ^ ® tr ^z_9tz^ — z^_^dtz^^ dt = — ^ lu tTz^_^z^ (22) 

Note that u > and v > 0. It is straightforward to 
show that ( [l9| ) equals ismX{u — v), whereas ( pO| ) equals 
(cos A — l)(u + t;), which completes the proof of (|]). 

Now we consider a single channel pump, S{t) e U{2). 
In this case, z and are complex numbers. For harmonic 
driving signal, without loss of generality, one can write 



z(t) = ziVi coa{nt + 9) + Z2V2 cos{nt), 



(23) 



where Vi^2 are pumping signal amplitudes, and complex 
parameters zi^2 depend on microscopic details. From ( |2l| ) 
we find the Poisson rates 



1 



1 



ziVie-'<^ + Z2V2\ (24) 



Note that u and v vanish at particular signal amplitudes 
ratio V1/V2 and phase 6. Once the two Poisson processes 
(^ are reduced to one, the current to noise ratio gives 
elementary charge, I/J = ±e~^. This happens at the 
extrema of I/J as a function of w = iVi/V2)e'^^ , for ( |2^ ) 
reached sd w — —Z2/ zi, —Z2/z\. This behavior is illus- 
trated in Fig.^ where Z2/ z\ = 0.4 -I- 0.8z is used. 

Reducing the counting statistics (^ to purely Poisso- 
nian by varying pumping signal parameters, in principle, 
is possible for any number of channels n. However, since 
the number of parameters to be tuned is 2n^, this method 
is practical perhaps in the single channel case only. Al- 
though the method is demonstrated for non-interacting 
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